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Abstract 

A closed expression for the differential cross section of the large-angle Bhabha 
e + e~ scattering which explicitly takes into account the leading and next-to-leading 
contributions due to the emission of two hard photons is presented. Both collinear 
and semi-collinear kinematical regions are considered. The results are illustrated by 
numerical calculations. 
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1 Introduction 



The large-angle Bhabha process is well suited for the determination of the luminosity C at 
e + e~ colliders of the intermediate energy range y/s = 2e ~ lGeV [1,2]. Small scattering 
angle kinematics of Bhabha scattering is used for high-energy colliders such as LEP I 
[3]. As far as 0.1% accuracy is desirable in the determination of C, the corresponding 
requirement 
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< i<r 3 (l) 



on the Bhabha cross section theoretical description appears. The quantity Aa is an un- 
known uncertainty in the cross section due to higher order radiative corrections. A great 
attention was paid to this process during the last decades (see review [4] and references 
therein). The Born cross section with weak interactions taken into account and the first 
order QED radiative corrections to it were studied in detail [5] . Both contributions, the one 
enhanced by the large logarithmic multiplier L = ln(s/m 2 ) (where s = (p + + p-) 2 = 4e 2 
is the total center-of-mass (CM) energy squared, m is the electron mass), and the one 
without L are to be kept in the limits (1): aL/n, a /it. As for the corrections in the second 
order of the perturbation theory, they are necessary in the leading and next-to-leading 
approximations and take the following orders, respectively: 



-)*L>, ( a -\\. (2) 



irj \irj 

The total two-loop (~ (a/n) 2 ) correction could be constructed from: 1) the two-loop 
corrections arising from the emission of two virtual photons; 2) the one-loop corrections 
to a single real (soft and hard) photon emission; 3) the ones arising from the emission of 
two real photons; 4) the virtual and real e + e~ pair production [6]. As for the corrections 
in the third order of perturbation theory, only the leading ones proportional to (aL/n) 3 
are to be taken into account. 

In this paper we consider the emission of two real hard photons: 

e + {p + ) + e"(p_) -> e+(g + ) + e"(g_) + 7(^1) + 7(^2)- (3) 

The relevant contribution to the experimental cross section has the following form 

a cxp = J d<7 + 0_, (4) 



where 0+ and 0_ are the experimental restrictions providing the simultaneous detection 
of both the scattered electron and positron. First, this means that their energy fractions 
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should be larger than a certain (small) quantity s^/e, e t h is the energy threshold of the 
detectors. The second condition restricts their angles with respect to the beam axes. They 
should be larger than a certain finite value ipo ( ipo ~ 35° in the experimental conditions 
accepted in [1]): 

7T - ip > 6L, 0+ > ijj , 0± = q^p_ , (5) 



where 6± are the polar angles of the scattered leptons with respect to the beam axes 
(?>_)• We accept the condition on the energy threshold of the charged-particle registration 
<Z± > £th- Both photons are assumed to be hard. Their minimal energy 



uj min = Ae, A < 1, (6) 
could be considered as the threshold of the photon registration. 

The main (~ (aL/n) 2 ) contribution to the total cross section (5) comes from the collinear 
region: when both the emitted photons move within narrow cones along the charged parti- 
cle momenta (they may go along the same particle). So we will distinguish 16 kinematical 
regions: 

aki and ak 2 < , aki and bk 2 < 9 , 

Tfl 

— «0 O «1, a^b, a,b = p_,p + ,q_,q + . (7) 

e 

The matrix element module square summed over spin states in the regions (7) is of 
the form of the Born matrix element multiplied by the so-called collinear factors. The 
contribution to the cross section of each region has also the form of 2 — > 2 Bhabha cross 
sections in the Born approximation multiplied by factors of the form 



da™ 11 = da 0i 



aKx.y^ln^^+b^x.y^ln^ 



HI 2 



(8) 



where Xj = Wj/e, y\ = q°_/e, y 2 = q+/s are the energy fractions of the photons and of 
the scattered electron and positron. The dependence on the auxiliary parameter # will 
be cancelled in the sum of the contributions of the collinear and semi-collinear regions. 
The last region corresponds to the kinematics, when only one photon is emitted inside the 
narrow cone Q\ < 6$ along one of the charged particle momenta. And the second photon 
is emitted outside any cone of that sort along charged particles (02 > ): 

d< = ^ln(^)d^(k 2 ), (9) 

where daj; has the known form of the single hard bremsstrahlung cross section in the 
Born approximation [7]. 
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Below we show explicitly that the result of the integration over the single hard photon 
emission in eq. (9) in the kinematical region 9\ > Oq (9\ is the emission angle of the second 
hard photon with respect to the direction of one of the four charged particles) has the 
following form 

J d^(k 2 ) = -21n(|)a i (x,y)d ( jj ) + da 1 . (10) 

The collinear factors in the double bremsstrahlung process were first considered in pa- 
pers of the CALCUL collaboration [8]. Unfortunately they have a rather complicated 
form which is less convenient for further analytical integration in comparison with the 
expressions given below. The method of calculation of the collinear factors may be con- 
sidered as a generalization of the quasi-real electron method [9] to the case of multiple 
bremsstrahlung. Another generalization is required for the calculations of the cross section 
of process e + e~ — > 2e + 2e~ [6]. 

It is interesting that the collinear factors for the kinematical region of the two hard 
photon emission along the projectile and the scattered electron are found the same as for 
the electron-proton scattering process considered by one of us (N.P.M.) in paper [10]. 

There are 40 Feynman diagrams of the tree type which describe the double bremsstrahlung 
process in e + e~ collisions. The differential cross section in terms of helicity amplitudes 
was computed about ten years ago [8,11]. It has a very complicated form. We note that 
the contribution from the kinematical region in which the angles (in the CM system) 
between any two final particles are large compared with m/e is of the order 

7T Z E Z 



(r is the classical electron radius). So, the corresponding events will possess poor statistics 
at the colliders with the luminosity £ ~ 10 31 — 10 32 cm~ 2 s _1 . More probable are the cases 
of double bremsstrahlung imitating the processes e + e~ — > e + e~ or e + e~ — > e + e~7, which 
corresponds to the emission of one or two photons along charged-particle momenta. 



2 Kinematics in the collinear region 



It is convenient to introduce, in the collinear region, new variables and transform the 
phase volume of the final state in the following way (from now on we will work in the CM 
system) : 

f i-p /- d 3 q_ d 3 q+d 3 k! d 3 k 2 (4) 
J J 16q_q+u;iCc>2(27r) 8 
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Table 1 



rji and Aj for different collinear kinematics. 







9-9- 


P+P+ 


9+9+ 


P-P+ 


9-9+ 


P-Q- 


P+9+ 


p-q+ 


P+Q- 


m 


y 


1 


1 


1 


1 — Xl 


1 


1 — Xl 


1 


1 — Xl 


1 


m 


1 


1 


y 


1 


1 — x 2 


1 


1 


1 — X\ 


1 


1 — Xl 


Ai 


1 


l 

y 


l 


1 


1 


l 

1 — Xl 


1 + aa 

2/1 


1 


1 


1 + ^2 
2/1 


A 2 


l 


1 


l 


1 

y 


1 


1 

1—X2 


1 


1+21 
2/2 


l + a 

2/2 


1 



4 2 1 1 27r 

WIT 



dxi / dx 2 xix 2 / — / dzi / dz 2 / dr q , (12) 



4(27r) 6 7 V V 2tt 

A A 
dFq = J 4q°q°(27rV (4)fap - + %P+ " AlQ " " A2q+) ' 



4q0q + (27r) 5 

^1,2 = , <P = k 1± _k 2 ^ Xi = — , 2; = — > 1, A = , 

\ m J e \ m J e 

where 6i (i = 1, 2) is the polar angle of the i-photon emission with respect to the mo- 
mentum of the charged particle that emits the photon; rj± and A± depend on the specific 
emission kinematics, they are given in Table 1. 



The columns of the Table correspond to a certain choice of the kinematics in the following 
way: p_p_ means the emission of both the photons along the projectile electron, 
means that the first of the photons goes along the projectile positron; the second, along 
the scattered electron, and so on. The contributions from 6 remaining kinematical regions 
(when the photons in the last 6 columns are interchanged) could be found by the simple 
substitution x\ <-> x 2 . We will use the momentum conservation law 

ViP- + V2P+ = Ai<7_ + A 2 g+ , (13) 



and the following relations coming from it: 



Vi+V2 = Aij/i + \ 2 y 2 , Ait/i sin 0_ = \ 2 y 2 sin 6>+, 



A 2 y 2 = 



vi 



rjl)c 



1/1,2 = — , 



V1+V2 + (rj2 -Vi)c 



(14) 



Each of 16 contributions to the cross section of process (3) can be expressed in terms of 
the corresponding Born-like cross section multiplied by its collinear factor: 
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2!V27r7 2 M) 

dcr„fa, A) = — A) dlfa, A), Bfa, A) = (^±£±- St "' " 
s V st 



dli(?7, A) = q Q d Q q+ ^ (4) fap, + ry 2 p + - Aig_ - A 2 g+) 



4irr)ir) 2 dc 



*i*![c(»72 -»7i) + ^] 2 ' 

20 z 2tt 

K( V ,X)=m 4 J dz! y dz 2 I ^/C^A), 







r / x x2 ?7i(l-c) 

t = {ViP- ~ MQ-) = 



V1+V2 + (r)2 ~ Vi) c ' 
5 = (?7iP- + W+) 2 = ^ 2 ViV2 = s^i^2, s + t + u = 0. 

The sum over (77, A) means the sum over 16 collinear kinematical regions. The correspon- 
ding (77, A) could be found in Table 1. The quantities /C^ (77, A) are as follows: 

}C(p_p_) = -A(A 1 , A 2 , A, Xl ,x 2 , y), /C(g_g_) = 2yA(B 1 , B 2 , B, -), 

y y y y 

2 £^ j. 1 

£(P+P+) = ~A(Ci, C 2 , C, x u x 2 , y), JC(q+q+) = 2yA(D u D 2 , D, -, -, -), 

y y y y 

A( A A A \ yA2 yAl 1+y2 r ? + W2 

A(A U A 2 , A, Xl , X2 ) = -— - — + + - n — 

, r^ + yr 1 2m 2 (y 2 + rj) 2m 2 (y 2 + r%) 



}C(p- P+ ) = 2K 1 K 2 , /C(p_g + ) = -2K X K Z , /C(p+g_) = -2K 4 ^ 5 , (17) 
/C(g_g + ) = 2K 6 ^ 7 , /C(p_g_) = -2K x K b , K,{p + q + ) = -2K,K 3 , 
gi 2m 2 g 2 2m 2 g± 2m 2 

1 = I^ + ^f' 2 = c 2 ~x^ r 2 + ~cf 1 z = D^t 2 ~~b\' 

gi 2m 2 #3 _ ^ Pi _ 

4 Cixin C 2 ' 5 B 2 x 2 t, B 2 ' 6 B lXl B\ ' 

g 2 2m 2 

Pi = l + r%, g 2 = l+rl, g 3 = yl + tl, g 4 = y\ + tj, 

ti = yi + x 2 , t 2 =y 2 + x 2 , y = l-X!-x 2 , 

yi, y 2 are the energy fractions of the scattered electron and positron defined in eq. (14). 
Expressions (17) agree with the results of paper [8] except for a simpler form of /C(g_g + ). 
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As for eq. (16) it has an evident advantage in comparison to the corresponding formulae 
given in paper [8]. Let us note that the remaining factors /C(p, q) could be obtained from 
the ones given in eq. (17) using relations of the following type: 



A. 


= (p- 


■ - fa) 2 


— m 2 , 
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(p_ -fa- k 2 ) 2 


— m 2 




= (?- 


+ fa) 2 


— m 2 , 


B = 


(g_ +fa + k 2 ) 2 


— m 2 , 




= (ki 




— m 2 , 
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(fa + fa -p+f 


— m 2 
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= (? + 


- + ki) 2 


— m 2 , 


D = 


(q+ + fa + fa) 2 


— m 2 



JC(p-q-)(x 1 ,x 2 ,A 1 ,B 2 ) = JC(q-P-)(x2,x u A 2 ,Bi). (18) 

Note also that terms of the kind m^j{B\C 2 ) do not give logarithmically enhanced con- 
tributions, and we will omit them below. The denominators of the propagators entering 
into eqs. (16), (17) are: 



(19) 



For further integration it is useful to rewrite the denominators in terms of the photon 
energy fractions x\^ 2 and their emission angles. In the case of the emission of both the 
photons along p_ we would have 



— - = — Xi(l + Zi) - X 2 (l + Z 2 ) + XiX 2 (zi + Z 2 ) + 2%iX 2 JZiZ 2 COS0, 

^ = -x t (l + Zl ), (20) 

where Zi = (e9i/m) 2 , <p is the azimuthal angle between the planes containing the space 
vector pairs (p_ , ki) and (p_ , k 2 ). In the same way one can obtain in the case fa , fa\\q- : 



-iL = ^1(1 + yf Zl ) + —(l + y 2 z 2 ) +x 1 x 2 (z 1 + z 2 ) + 2x 1 x 2v /ili^cos0, 
m 2 yi y 1 

4 = % + !/?*). (21) 
m 2 y 1 

Then we perform the elementary azimuthal angle integration and the integration over 
z\ , z 2 within the logarithmical accuracy using the procedure suggested in paper 



z z 2i 

a = m 4 J dzi J dz 2 J ^a. (22) 



The list of the relevant integrals is given in Appendix A. In this way one obtains the 
differential cross section in the collinear region: 
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a 4 L d 3 q + d 3 q_ dx l dx 2 , „ J 1 

do-con = — o — o 1 1 + ^1,2) S — 2 

47r 2 s q^q 1 ! x x x 2 [yri 

2 

+ (y 2 + ^1) In + ^1^2(2/ - a:ix 2 ) - 2ri5f 5 



■(L + 21)g lg5 



[Bp_p_5p_p_ + Bp + p + 5. 



(23) 



, - (L + 2/ + 4 In y)#i# 5 + (y 2 + r\) In + x x x 2 (y - x x x 2 ) - 2r 1 g 1 

yr( L2 x 2 y 



X [Bq-q-fiq-q- + -Sg + g + ^<j + g+] + Bp_ p+ 5p_p + 



(L + 2/)^-2^-2^ 

ri r 2 ri r2 



(L + 2Z + 21n(r 1 r 2 ))^-2^-2^ 
r\r 2 ri r 2 



[-^P— <?- ^P-q- Bp + q_5p + q_] 
~\~[Bp + q + 8p +q+ + -B p _q + (5p_g + ] 



(L + 2l + 21n yi )-^-2^-2-^ 

n 2/1*1 

(L + 2/ + 21n 2/2 )^^-2^-2-^ 4 



^12/2^2 n 1/2*2 



We used the symbol Pi ;2 for the interchange operator (Vi, 2 f(xi,x 2 ) = f(x 2 ,Xi) ). We 
used the notation (see also eq. (17)): 



I = In 



01 



2 1 2 
#5 = 2/ +r 1 , 



(24) 



where #0 is the collinear parameter. Delta-function S Ptq corresponds to the specific con- 
servation law of the kinematical situation defined by the pair p, q (see Table 1): 5 Ptq = 
S^(r] 2 p + +r]ip_ — \\q- — \ 2 q+). Besides, we imply that the first photon is emitted along 
the momentum p; and the second, along the momentum q (p, q = p_, p + , q_, q + ). These 
5-functions could be taken into account in the integration as is made in the expression 
for dl(?7, A) (see eq. (15)). Finally, we define 



VAit r] 2 s 



t=(p--q-) 2 . 



(25) 



3 Contribution of the semi— collinear region 



We will suggest for definiteness that the photon with momentum k 2 moves inside a narrow 
cone along the momentum direction of one of the charged particles, while the other photon 
moves in any direction outside that cone along any charged particle. This choice allows 
us to omit the statistical factor 1/2!. The quasireal electron method [9] may be used to 
obtain the cross section: 



da sc = 



or 



d 3 q_ d 3 q+ d 3 ki d 3 k 2 f JC P _ 



32S7T 4 



q°_q+K 



V 



k° 2 \p.k 2 5p - Rp - 
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We omitted the terms of the kind m 2 /(p-k 2 ) 2 in eq. (26) because their contribution does 
not contain the large logarithm L. The quantities entering into eq. (26) are given by: 



V = - : + 



s 



t' t 



kip + ■ hp- hq + ■ hq_ kip + ■ hq+ hp- ■ hq. 



u' u 



V is the known accompanying radiation factor; /Q are the single photon emission collinear 
factors: 



IC = IC K, = y± + ^— + ^ K, = V i + ^— + X ^ (28) 

P+ x 2 r 2 q ~ x 2 (yi + x 2 ) ' q+ x 2 (y 2 + x 2 ) 

Quantities Ri read: 



R p _ = R[sr 2 , tr 2 , ur 2 , s', t', u'], R p+ = R[sr 2 , t, u, s', t'r 2 , u'r 2 ], 

R q =R[s,t-,u,s'-,t',u'-l R q =R[s,t,u-,s'-,t'-,u'], (29) 

yi yi yi ' 2/2 2/2 1/2 

where the function R has the form [12]: 



R[s, t, u, s', t', u] = —^[ss\s 2 + s' 2 ) + tt'(t 2 + t' 2 ) + uu'(u 2 + u' 2 )}, 
s s tt 

S=(P++P-) 2 , s'=(q+ + q^) 2 , t=(p_-g_) 2 , 

t' = (p + -q + ) 2 , u=(p_-q + ) 2 , u' = (p + -q_) 2 . (30) 
Finally, we define 



5 p _=5^\p_r 2 +p + - q+ - q _-k 1 ), 
$p+ =5 (4) (P- +P+r 2 -q+-q--h), 
5 q _=^\p_+p + -q + -q_ V -±±^-k 1 ), 

S q+ = 5 {4) (p- + p+ - q+ V2 + X2 -q-~ h). 

1/2 



(31) 



Performing the integration over angular variables of the collinear photon we obtain 
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dcr s 



a 



l L d 3 q_ d 3 q + d 3 !^ 



16S7T 3 
1 



q°_q%l4 



dx 2 V /C p _[R p _<5 p _+R p+ £ 



p+j 



H JC q+ R q+ 5 q+ H }C q _R q _S q 

V2 V\ 



(32) 



To see that the sum of cross sections (23) and (32) 



/f](T SC 
dOi(^-) (33) 

does not depend on the auxiliary parameter 9 . We verify that terms L ■ I from eq. (23) 
cancel out with the terms 



f^n-L-l, (34) 
2tt J hqi y ' 



which arise from 16 regions in the semi-collinear kinematics. 



4 Numerical results and discussion 



We separated the contribution of the collinear and semi-collinear regions using the auxi- 
liary parameter #o. By direct numerical integration according to the presented formulae 
we had convinced ourselves that the total result is independent on the choice of 9 . 

It is convenient to compare the cross section of double hard photon emission with the 
Born cross section 



2 COS V>0 , 2X 2 

a^ = ^ [ 3 -±^] dc. (35) 



2s J V 1 - c 

- cos ipo 



For illustrations we integrated over some typical experimental angular acceptance and 
chose the following values of the parameters: 



rj} = vr/4, = 0.9 GeV, A 1 = 0.4, A = 0.05 , 6 = 0.05, 

L = 15.0, / = -7.38, (36) 

where Ai defines the energy threshold for the registration of the final electron and 
positron: q± > £ th = eAi. Note that restrictions on 9 (7) and (12) (z = exp{L + /} ^> 1) 
are fulfilled. 
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For the chosen parameters we get 



coll 

cr Born = 1.2 mkb, -5— • 100% = -0.25 % 

7 ^Born 

^sc ^.sc _|_ ^coll 



100% = 0.81 %, 5a tot = ' • 100% = 0.56 %. (37) 



0-Born ' ' 0-Born 



The phenomenon of negative contribution to the cross section from the collinear kine- 
matics is an artifact of our approach. Namely, we systematically omitted positive terms 
without large logarithms, among them we dropped terms proportional to I 2 . The can- 
cellation of I 2 terms can be seen only after adding the contribution of the non-collinear 
kinematics (when both photons are emitted outside narrow cones along charged-particle 
momenta). The non-collinear kinematics does not provide any large logarithm L. 

Both quantities ct co11 and a sc depend on auxiliary parameter 9 . We eliminated by hands 
from eq. (23) the terms proportional to I and obtained the following quantity: 



cr, 



bare 
coll 



^jBorn 



•100% = 1.43%. (38) 



This quantity corresponds to an approximation for the correction under consideration in 
which one considers only the collinear regions and takes into account only terms propor- 
tional to L 2 and L (all terms dependent on 9q are to be omitted). Having in mind the 
cancellation of #o~dependence in the sum of the collinear and semi-collinear contribu- 
tions, we may subtract from the value of the semi-collinear contribution the part which 
is associated with /: 

_bare 

= Cr sc + (ct co11 - <£™), -f^ ■ 100% = -0.87 %. 

Looking at bare quantities one can get an idea of relative impact of two considered regions. 
We see that at the precision level of 0.1% the next-to-leading contributions of semi- 
collinear regions are important. 

In figure 1 we illustrated the dependence on parameter A of the bare collinear contribution 
for different fixed values of Ai. Large growing in the region of small A corresponds to an 
infrared singularity, which will be cancelled after adding contributions of virtual and soft 
photon emission. 
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Appendix A 



We present here the list of integrals (see eqs. (19 - 22)): 



A, 



A 2 A X 



X\X 2 r\ 

u 

AA\ x\x 2 r x 



1 



1 x 2 r\ 
-Lq + In 

2 xiy 

1 x 2 r\ 
-L + In 

2 x x y 



1 + 



XiX 2 

y 



AA{ 



x\x 2 r\ ' 



(39) 



J 2 



1 



A) 



A X A 2 x\x 2 AtB 2 
L = In z = L + /, 



yix±x 2 



ff2 



-(L + 2\n yi ), 

As 2 
m z 



The remaining integrals could be obtained by simple substitutions defined in eqs. (19 - 
22). 
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